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INTRODUCTION 


Lidar  has  become  a  very  popular  instrument  for  use  in  the  study  of  atmos¬ 
pheric  extinction.  Lidar  theoiy,  however,  suffers  from  the  limitation  that,  while  there 
is  only  one  equation,  there  eu-e  two  unknowns.  Most  people  using  lidars  avoid  this 
problem  by  making  assumptions  about  the  relationship  between  these  unknowns,  or 
make  assumptions  about  the  cheu’acteristics  of  the  atmospheric  aerosols. 

One  assumption  that  is  sometimes  made  is  that  the  atmosphere  is  horizontally 
homogeneous.  A  m  're  frequently  used  {assumption  is  that  the  relationship  between 
backscatter  {uid  extinction  is  a  constant.  These  assumptions  are  often  not  justified. 

Paulson  (1986)  reports  that  lidar  returns  frequently  show  pronounced  horizon¬ 
tal  atmospheric  inhomogeneities.  Figure  1  shows  two  examples  of  this.  The  curves 
are  the  natured  logarithm  of  the  range-compensated  power  received  plotted  as  a  func¬ 
tion  of  range.  If  the  atmosphere  were  homogeneous  along  the  propagation  path,  these 
curves  would  be  straight  lines,  with  the  magnitude  of  the  slope  equal  to  twice  the 
atmospheric  extinction  coefficient. 

Workers  who  assume  constant  backscatter-to-extinction  ratios  usually  make  a 
guess  as  to  what  this  ratio  is  and  proceed  from  there.  Fenn  (1966),  however,  states 
that  theoretical  considerations  show  that  no  generally  valid  relation  exists  between 
the  atmospheric  extinction  coefficient  and  the  atmospheric  backscatter  properties.  He 
further  says  that  it  is  therefore  clear  that  such  a  relation  cannot  be  unique  and  will 
change  from  one  situation  to  another. 

Mulders  (1984)  states  that  the  relationship  between  backscatter  emd  extinc¬ 
tion  may  be  very  sensitive  to  the  chemical  composition  of  the  aerosol.  His  measure¬ 
ments,  he  says,  show  that  it  may  vary  during  one  day  in  one  place.  When  data  for  one 
month  were  considered,  no  relationship  between  backscatter  coefficient  and  extinc¬ 
tion  coefficient  was  evident. 

A  dual-lidar  method  for  measuring  extinction  coefficient  profiles  was  proposed 
by  Paulson  and  Powers  in  late  1985.  This  method  uses  two  lidars  set  at  opposite  ends 
of  a  propagation  path  and  pointed  at  each  other.  No  assumptions  need  to  be  made 
about  the  atmosphere  since  there  are  now  two  equations  available,  {md  the  extinction 
coefficient  profile,  integrated  extinction,  and  the  backscatter  coefficient  profile  can  be 
obtained.  The  theory  is  discussed  in  detail  by  Paulson  (1987)  and  by  Hu^es  and 
Paulson  (1988). 

Recently,  Kunz  (1987)  has  suggested  the  same  technique  with  a  somewhat  dif¬ 
ferent  mathematical  approach,  but  it  appears  that  he  has  not  tested  it  yet. 

The  purpose  of  this  report  is  to  describe  and  demonstrate  an  extension  of  this 
technique  to  obtain  optical  depth  measurements  in  the  vertic£il  direction.  While  the 
mathematical  derivations  have  been  documented  in  previous  reports,  they  are 
repeated  here  for  completeness. 
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Figurs  1.  S(R)  data  plotted  as  a  function  of  range  for  Ildar  #025091  taken  on  27  August  1985. 
Visibility  was  about  10  km. 
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MATHEMATICAL  DERIVATIONS 


EXTINCTION  COEFFICIENT  PROFILES 


The  quantity  S(R)  is  the  natural  logarithm  of  the  backscattered  power 
received  from  range  R,  multiplied  by  R^.  If  the  two  lidars  are  separated  a  distance  d, 
as  showm  in  figure  2,  and  the  origin  is  at  lidar  I,  the  equation  for  S(R)  for  the  first 
lidar  is 


S.W  =/«(Cn)+/«[^W]-  air)dr 

Jo 


and  that  for  the  second  lidar  is 


(1) 


52(/?)=M(Ca)+;/i(y?(/?)l-2  a(r)dr  (2) 

.1^ 

where  C,,  eind  C,2  are  the  instrumentation  constants  for  each  of  the  lidars,  a{rj  is  the 
extinction  coefficient  at  range  r,  and  fi(R)  is  the  backscatter  coefficient  at  range  R. 


If  equation  (2)  is  subtracted  from  equation  (1)  we  get 


S,(/?)-52(/?)  =/n(Cn-/n(C, 


-r 


a(r)dr  + 


2  [  a{r)dr  . 
J/t 


Since 


(3) 


[  a(r)dr  =  f  a{r)dr  -  f  o(r)dr  , 
'r  Jo  Jo 


equation  (3)  becomes 


(4) 


Si(R)-S2W  =ln(Cu-ln(Cn)-4  o(r)dr  +  2j^  a(r)dr  . 


(5) 


Taking  the  difference  between  equations  (1)  and  (2)  eliminates  the  backscatter 
coefficient.  Taking  the  derivative  of  equation  (5)  eliminates  the  requirement  that  con- 
st£mts  Cii  and  Ci2  be  known,  and  we  get 


dSi(R) -dSiiR)  =  -  4a(R)dR 


(6) 


or 


a(R)  = 


dSzCR)  dSdH) 
dR  dR 
4 


(7) 


The  calibration  curves  for  each  of  the  lidar  receivers  are  still  needed,  however, 
and  must  be  accurately  known  since  they  affect  the  slope  chEU’acteristics  of  the  S(R) 
curves.  Since  the  propagation  for  the  two  lidars  is  in  opposite  directions  with  respect 
to  the  origin,  the  slopes  of  S^(R)  and  SzfR)  should  have  opposite  signs  under  homoge¬ 
neous  conditions. 
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backscatter  coefficient  profiies. 


INTEGR.^TED  EXTINCTION  OR  OPTICAL  DEPTH 


While  equation  (7)  will  give  good  extinction-coefficient  profiles  under  many 
conditions,  small-scale  irregularities  in  the  atmospheric  aerosols  may  be  seen  some¬ 
what  differently  by  the  two  lidars.  This  can  cause  brief  positive  excursions  in  the 
S(R)  difference  curve,  which  can  appear  as  negative  extinction  coefficients.  In  a 
ground-based  system  this  condition  can  be  minimized  by  careful  alignment  of  the  two 
lidars  so  that  they  su'e  looking  along  exactly  the  same  path;  but  in  the  case  of  a  verti¬ 
cal  system,  where  one  lidar  is  in  an  aircraft,  this  is  not  so  easy.  It  appears,  however, 
that  the  integrated  extinction,  or  optical  depth,  obtained  from  the  S(R)  difference 
curve  is  not  so  sensitive  to  these  positive  excursions. 

If  equation  (7)  is  integrated  between  two  distances,  i?,  and  i?2,  we  get 


a(r)dr  = 


4 


(8) 


Therefore,  just  taking  the  difference  in  the  S(R)  difference  curve  at  Range(l)  gmd 
Range(2)  and  dividing  by  4  gives  the  integrated  extinction  over  that  distance. 


LIDAR  CALIBRATION 

The  lidars  used  in  these  measurements  were  Visioceilometer  lidars  operating 
at  1.06  (i,m  wavelength. 

Lidar  2  (#025091)  was  calibrated  with  the  assistance  of  personnel  at  TNO  in 
the  Netherlands  (Ferguson  &  Paulson,  1986).  It  was  subsequently  recalibrated  along 
with  lidar  1  (#025090)  using  neutral  density  filters.  Neutral  density  filters  from  0  to 
70  dB  in  5-dB  steps  were  placed  in  front  of  the  lidar  receiver.  Several  shots  were 
made  at  each  attenuation  and  these  were  averaged.  A  fifth-order  fit  was  made  for  the 
recorded  signal  eunplitude  versus  the  attenuation  for  each  of  the  lidars. 


LIDAR  COMPARISON 

The  two  lidars  were  set  up  side  by  side  overlooking  the  oceeui  at  an  elevation 
of  about  30  m.  Elevation  an^e  for  each  was  adjusted  so  that  the  cross  hair  rested  on 
the  horizon.  The  vertical  cross  hair  for  each  lidar  was  aligned  on  the  mast  of  a  sail¬ 
boat  about  4  or  5  miles  away  to  make  sure  that  they  were  edigned  in  azimuth.  The 
two  lidars  were  then  fired  nearly  simultaneously  for  a  number  of  shots.  Figure  3 
shows  a  compeu'ison  of  the  S(R)  curves  for  two  of  these  shots.  These  show  quite  good 
agreement  out  to  about  1  km.  Background  n^'ise  becomes  important  beyond  that. 


PftTA  SET  5  02/17/89  Lidar  IJ025091(Red)  #025090 (Green) 


0  .5  1 

RANGE  (kn) 


DATA  SET  6  02/17/89  Lidar  «025091(Red)  #025090 (Green) 


0  .5  1 


RANGE  (kn) 

Figure  3.  Comparison  of  the  two  lidars  f<red  along  a  parallel  horizontal  path  over 
the  ocean.  The  increasing  S(R)  with  range  would  indicate  a  gradient  in  the 
aerosol  distribution. 
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VERTICAL  MEASUREMENTS 

EQUIPMENT  SETUP 

The  first  lideir  was  mounted  on  a  tripod  and  pointed  in  the  vertical  direction. 
The  location  of  this  lidfir  was  at  NOSC  building  15  at  the  south  end  of  Point  Loma. 
This  site  was  about  25  to  30  m  above  sea  level.  The  second  lidar  was  installed  in  a 
Piper  Nav^o  aircraft  and  pointed  vertically  down  through  a  hole  in  the  floor. 

The  aircraft  made  repeated  passes  over  the  ground  site  at  a  constant  altitude 
of  2500  ft,  or  about  760  m  above  the  ground  site.  Personnel  on  the  ground  directed 
the  aircraft  via  radio  and,  when  it  was  directly  over  the  ground  site,  called  “mark,” 
and  both  lidars  were  fired.  The  data  were  recorded  on  Memodyne  digital  tape  record¬ 
ers  for  later  processing. 

DATA  ANALYSIS  AND  RESULTS 

The  first  set  of  data  was  taken  on  30  March  1989.  Ten  passes  were  made  at 
about  3-  or  4-min  intervals,  resulting  in  eight  data  samples.  The  lidar  in  the  aircraft 
failed  to  fire  the  first  two  passes. 

The  upper  graph  in  figure  4  shows  an  example  of  the  range-compensated 
returns  for  the  two  lidars.  The  lower  graph  shows  the  same  thing,  but  with  a  seven- 
point  running  average.  Minimum  altitude  for  the  ground  lidar  was  about  100  m.  This 
is  the  point  where  the  laser  beam  and  the  receiver  field-of-view  overlap.  The  near 
point  for  the  aircraft  lidar  was  about  680  m,  or  about  90  m  below  the  aircraft. 

Figure  5  shows  the  s£ime  thing  for  the  following  pass.  Data  Set  4.  In  this  case 
atmospheric  conditions  at  the  ground  changed  in  such  a  way  that  the  digitizer  lim¬ 
ited  for  return  signals  closer  than  about  120  m.  (This  point  is  marked  with  an  arrow 
in  the  upper  graph.)  This  also  had  the  effect  of  reducing  signal  strength  fi’om  higher 
altitudes  so  that  signal  returned  from  greater  than  450  m  altitude  was  in  the  back¬ 
ground  noise.  This  same  limiting  occurred  for  Data  Sets  5  and  6. 

S(R)  difference  curves  8u*e  shown  in  figure  6  for  the  seven-point  running  aver¬ 
ages  for  these  two  data  sets.  Theoretically  the  slope  of  these  curves  should  never  go 
positive,  which  would  imply  negative  extinction.  While  the  slope  for  Data  Set  3  is 
generally  negative,  there  are  some  positive  perturbations  on  it.  These  are  the  result 
of  very  small  irregulju'ities,  which  are  seen  by  the  two  lideu^  somewhat  differently. 
Averaging  reduces  this  effect  somewhat,  but  cannot  completely  eliminate  it.  In  the 
S(R)  difference  curve  for  Data  Set  4,  the  curve  goes  horizontal  at  about  450  m  and 
then  starts  to  increeise.  This  is  the  point  where  the  ground  lidar  return  goes  into  the 
noise.  Below  about  120  m  the  digitizer  limited  for  the  ground  lidar,  otherwise  the  dif¬ 
ference  curve  might  have  been  somewhat  higher  in  this  region. 

If  we  apply  equation  (8)  to  the  S(R)  difference  curve  for  Data  Set  3,  we  get  an 
optical  depth  of  0.20  between  the  altitudes  of  120  euid  690  m.  Doing  the  same  thing 
ior  Data  ^t  4,  we  get  an  optical  depth  of  0.44  between  130  and  460  m.  This  has  been 
done  for  Data  Sets  5  through  9  as  well  and  the  results  shown  in  table  1.  Graphs  of 
these  data  are  in  the  appendix. 
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ALTITUDE 


FigurQ  4.  Upper  graph  is  range-compensated  lidar  returns  as  a  function  of  altitude  for 
the  two  lidars  for  Data  Set  3.  Lower  graph  is  seven-point  running  averages  for  the 
data  shown  in  the  upper  graph. 


DATA  SET  4  03/^0/^3 


Ground 


Aircraft 


ALTITUDE 

(km) 

DATA  SET  4 

03/30/^9 

Ground 


Aircraft 


ALTITUDE  (kn) 

Figure  5.  Upper  graph  is  range-compensated  lidar  returns  as  a  function  of  altitude  for  the 
two  lidars  for  Data  Set  4.  The  arrow  shows  the  altitude  below  which  the  digitizer  limited. 
Lower  graph  is  seven-  point  running  averages  for  the  data  shown  in  the  upper  graph. 
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DATA  SET  3  03/30/89 
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Figura  8.  S(R)  difference  curves  for  Data  Sets  3  and  4.  In 
return  from  the  ground-based  lidar  goes  into  the  noise  at 


Table  1.  Integrated  extinction  between  the  altitudes  indicated  for  each  of  the  data  sets.  Outside 
of  these  altitudes  the  signal  of  one  of  the  lidars  was  either  in  the  noise  or  the  digitizer  clipped. 

DATA  SET  ALTITUDE(m)  to  ALTITUDE(m)  OPTICAL  DEPTH 


3 

120 

690 

0.20 

4 

130 

460 

0.44 

5 

135 

590 

0.30 

6 

180 

465 

0.26 

7 

160 

520 

0.16 

8 

160 

460 

0.32 

9 

160 

550 

0.38 

In  Data  Sets  5  and  9  the  return  from  the  airborne  lidar  dipped  down  to  th.j 
noise  level  near  300  m  altitude.  This  caused  a  large  positive  excursion  in  the  S(R)  dif¬ 
ference  curve  at  that  point.  This  probably  does  not  ^ect  the  value  of  the  opticsil 
depths  calculated  since  the  signal  is  above  the  noise  at  the  points  used  for  the 
calculation. 

DISCUSSION  OF  RESULTS 

The  lidars  used  have  about  a  70-dB  dynamic  range.  While  this  appears  to  be 
quite  large,  variations  in  atmospheric  conditions  can  greatly  exceed  it,  especially 
when  making  measurements  in  the  vertical.  In  the  case  of  the  ground-ba^  lidar, 
hi^  concentrations  of  aerosols  close  to  the  ground  cause  strong  signal  return  from 
close  to  the  lidar,  sometimes  exceeding  the  limits  of  the  digitizer.  These  same  hi^ 
concentrations  tend  to  reduce  signal  returned  from  higher  up,  where  there  are  fewer 
aerosols  and  the  backscattered  signal  is  weaker.  This  condition  works  in  favor  of  the 
airborne  lidar.  Since  the  aircraft  is  at  altitudes  where  the  aerosol  concentrations  are 
low,  returns  from  these  aerosols  are  close  to  the  lidar,  and  the  signal-to-noise  ratio  is 
good.  At  lower  altitudes,  farther  from  the  lidar,  the  aerosol  concentration  increases, 
causing  an  increase  in  the  strength  of  the  backscattered  signal  and  thus  maintaining 
a  good  signal-to-noise  ratio  in  most  cases.  The  reverse  may  also  occur  when  there  6tre 
stratus  conditions.  Then  the  airborne  lidar  is  in,  or  near,  hi^  aerosol  concentrations 
and  the  one  on  the  ground  may  not  be. 

AIRCRAFT  CLIMBING  SPIRAL  TESTS 

On  completion  of  the  dued-lidar  measurements,  the  aircraft  was  flown  to  a 
location  near  the  Coronodo  Islands.  There  a  climbing  spiral  was  made  from  100  to 
5000  ft  while  instruments  measured  dew  point,  temperature,  and  aerosol  size  distri¬ 
butions.  At  1000-ft  intervals  the  aircraft  leveled  momentarily,  and  the  lidar  was  fired 
vertically  downward. 
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LIDAR  S(R)  VERSUS  RELATFSTE  HUMIDITY 

The  lidar  shots  were  combined  to  generate  an  S(R)  profile  as  a  function  of  alti¬ 
tude  up  to  the  5000-ft  level.  This  curve  is  shown  in  figure  7.  Dew  point  and  tempera¬ 
ture  were  used  to  calculate  relative  humidity  as  a  function  of  altitude,  and  it  and  the 
S(R)  profile  were  compared.  This  comparison  is  shown  in  figure  8.  The  two  curves 
follow  each  other  quite  well  from  the  surface  up  to  about  380  m.  Above  this  altitude 
the  relative  humidity  remains  fairly  constant  at  around  20  percent,  but  the  S(R)  goes 
to  a  minimum  of  -8.1  then  increases  again  to  -6.2  before  gradually  decreasing  again. 
There  does  appear  to  be  a  tendency  for  the  two  curves  to  follow  each  other  up  to  at 
least  800  meters,  however.  There  also  is  the  possibility  of  aerosols  at  these  altitudes 
which  are  not  affected  by  relative  humidity. 

LIDAR  S(R)  VERSUS  KNOLLENBERG  S(R) 

Aerosol  sizes  were  measured  continuously  from  100  to  5000  ft  using  a  Knol- 
lenberg  Instrument.  These  measurements  were  taken  in  4-s  sample  intervals,  and  size 
distributions  were  calculated  for  each.  A  Mie  program  was  used  to  calculate  backscat- 
ter  and  extinction  coefficients  for  each  of  these  samples.  The  resulting  backscatter 
and  extinction  coefficients  were  then  used  in  a  lidar  equation  to  calculate  an  S(R) 
curve  as  a  function  of  altitude.  This  curve  is  shown  in  figure  9  and  compared  to  that 
measured  with  the  lidar.  The  two  curves  show  very  good  agreement  up  through  the 
inversion.  Above  that,  the  lidar  S(R)'s  are  somewhat  higher  than  those  calculated 
from  the  Knollenberg  measurements.  Also,  the  Knollenberg  values  appear  to  fluctuate 
much  more.  This  is  because  the  Knollenberg  counts  very  few  aerosols  in  a  sampling 
period  at  these  altitudes,  so  the  measured  distribution  is  probably  not  a  good  indica¬ 
tion  of  the  true  size  distribution  for  these  conditions. 

CONCLUSIONS 

Under  some  conditions  the  dual-lidar  method  can  be  used  to  get  measured 
optical  depths  in  the  vertical  without  making  assumptions  about  the  chmacteristics 
of  the  atmospheric  aerosol  distributions.  For  this  to  be  done  successfully,  both  lidars 
must  have  good  signal-to-noise  ratios  over  the  height  remge  of  interest.  At  the  same 
time  the  return  signal  cannot  be  so  strong  as  to  limit,  or  have  the  digitizer  clip,  at 
least  in  the  range  of  interest.  Also,  some  data  smoothing  is  required  to  reduce  the 
effects  of  small  irregularities,  which  may  be  seen  somewhat  differently  by  the  two 
lidars. 

Range-compensated  lidar  returns  and  relative  humidity  variations  showed  the 
same  trends  up  through  the  inversion.  Even  above  the  inversion  there  was  some  cor¬ 
respondence.  However,  it  appeeu’s  that  there  may  be  other  aerosols  contributing  to 
the  lidar  return  in  this  region  that  are  not  dependent  on  relative  humidity. 

The  range-compensated  lidar  returns  showed  good  agreement  with  an  S(RJ 
curve  calculated  with  the  use  of  Knollenberg  aerosol  size  distributions  up  through  the 
inversion.  It  appears,  however,  that  above  the  inversion  the  number  of  aerosols  inter¬ 
cepted  by  the  Knollenberg  weis  not  sufficient  to  provide  a  good  measure  of  the  true 
size  distributions. 
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Figuro  8.  Comparison  of  range-compensated  lidar  return  as  a  function  of  altitude  with  relative 
humidity  plotted  as  a  function  of  altitude. 
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APPENDIX 

S(R)  CURVES  AND  S{R)  DIFFERENCE  CURVES 
FOR  DATA  SETS  6  THROUGH  9 
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Figurs  A-2.  S(R}  curves  for  for  Data  Set  6.  Lower  curves  are  seven-point 
running  averages  of  the  upper  curves. 
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Figur«  A-5.  S(R)  curves  for  for  Data  Set  9.  Lower  curves  are  seven-point 
running  averages  of  the  upper  curves. 
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Figure  A-8.  S(R)  difference  curve  for  Data  Set  9. 
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